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NOMENCLATURE 


A  the  blocked  area  between  the  body  and  duct  wall 

A  area  of  duct 

o 

Cp  pressure  coefficient 

ACp  pressure  coefficient  correction 

ACp^  components  of  ACp  (i  =  1,  2,  ...) 

Cf  body  skin  friction  coefficient 

b 

Cq  suction  coefficient 

D  maximum  diameter  of  the  body 

max 

I  the  momentum  integral,  Eq.'s  (9)  and  (10) 

L  total  body  length 

£  length  from  forward  stagnation  point  of  body  to  beginning  of  duct 

SLq  distance  from  nose  of  body  to  point  on  duct  wall  where  turbulent 

boundary  layer  begins 

r  radius  coordinate 

free-stream  velocity  in  duct  far  upstream  of  body 

u  local  fluid  velocity 

uq  reference  velocity  in  potential  flow  computation 

x  axial  coordinate  where  origin  Js  upstream  from  nose  of  the  body 

x  body  axial  coordiante  (x  - 

x'  dummy  variable  of  integration 

6  boundary- layer  thickness 

v  kinematic  viscosity 

p  fluid  mass  density 

T  wall  shear  stress 

independent  variable  in  Equations  (15a)  and  (15b) 


e 
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Special  Subscripts 
b  refers  to  body 

1  refers  to  inviscid 

V  refers  to  viscous 

w  refers  to  duct  wall 
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INTRODUCTION 


The  pressure  distribution  on  a  body  of  revolution  will  be  altered 
from  that  observed  under  free  stream  conditions  when  the  body  is  operated 
in  a  closed  conduit  such  as  the  test  section  of  a  wind  or  water  tunnel. 

The  differences  are  a  result  of  blockage  effects  due  to  the  walls  of  the 
test  section  and  to  viscous  effects  originating  on  both  the  walls  and 
the  body  under  investigation.  As  pointed  out  by  this  author  [1,2]*,  the 
solid  blockage  effects  can  be  accurately  predicted  with  currently  used 
potential  flow  computer  codes  [3  through  7],  while  the  viscous  effects 
must  be  accounted  for  through  subsequent  analysis  of  the  potential  flow 
results.  Such  analysis  was  presented  in  References  [1]  and  [2]  under 
the  assumption  that  the  boundary  layer  flow  evolved  in  a  natural  flat 
plate-like  manner  over  the  body  and  duct  walls. 

In  some  instances,  the  body  under  investigation  may  have  a  form 
of  boundary  layer  control  such  as  heating,  cooling,  suction,  or  blowing. 

Under  these  conditions,  the  boundary  layer  growth  and  the  distribution 
of  wall  shear  stress  will  be  significantly  different  from  those 
calculated  or  measured  without  boundary  layer  control.  An  extension 
of  the  original  analysis  [1,2]  is  therefore  required  in  order  to  allow 
for  the  modified  boundary  layer  development  over  bodies  that  use  boundary 
layer  control.  This  report  gives  the  details  of  this  extension.  It  is 
assumed  that  the  boundary  layer  velocity  profile  and  wall  shear  stress 
distributions  are  available  as  input  for  this  analysis.  Hoffman  [8] 
has  developed  general  computational  codes,  based  on  the  work  of  Cebeci 
and  Smith  [9],  to  predict  these  parameters  for  bodies  with  roughness, 
cooling,  heating,  suction,  or  blowing.  As  a  numerical  example,  the 
analysis  presented  herein  will  be  used  to  predict  the  in-duct  body  pressure 
distribution  of  a  typical  axisymmetric  body  with  and  without  surface  suction. 


ANALYSIS 


Given  the  body  shape  and  cylindrical  duct  diameter  (2r  ),  it  is 
relatively  straightforward  to  calculate  the  inviscid,  in-duct  body 
pressure  distribution  [3  througti  7].  This  potential  flow  pressure 
distribution  is  denoted  by  Cpj-(x),  where  x  is  the  body  axial  coordinate, 
Figure  (la).  As  discussed  in  the  previous  papers  [1,2],  Cp_  should  be 
corrected  to  the  reference  velocity,  U^  which  would  occur  at  the  beginning 
of  an  assumed  infinitely  long  duct.  For  finite  distances  (£)  from  the 
nose  of  the  body  upstream,  the  flow  velocity  is  u  which  is  slightly 
higher  than  U^.  For  most  potential  flow  calculations  of  this  type, 

29,  +  L  represents  the  duct  length,  where  L  is  the  body  length  and 
then  „ 


cPl 


1  -  (i  -  Cp)  & 

•  ocv 


(1) 


In  Equation  (1),  Cp(x)  are  the  pressure  coefficient  values  computed 
by  the  numerical  code  being  used. 


» 


*Numbers  in  brackets  refer  to  references  at  end  of  report. 
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For  the  case  of  viscous  flow  in  the  duct,  boundary  layers  will 
develop  on  the  body  and  on  the  duct  wa_lls.  The  origin  of  duct  boundary 
layer  flow  is  assumed  known  and  is  at  x  =  as  schematically  expressed 
_in  Figure  (lb).  Note  also  that  a  new  coordinate  (x)  originates  at 
x  =  -A  and  that  the  velocity  at  x  =  0  is  Um.  A  goal  of  the  present 
analysis  is  to  derive  expressions  for  a  pressure  coefficient  correction 
(ACp)  such  that  Cpj  may  be  conveniently  corrected  for  the  viscous  flow 
effects,  i.e., 


CpV  =  Cpi  +  ACp 


(2) 


This  derivation  was  carried  out  in  References  [1]  and  [2]  for  assumed 
flat-plate  like  boundary  layer  development,  where  now,  the  analysis 
will  be  modified  to  account  for  the  non-conventional  types  of  boundary 
layer  development  that  occur  on  bodies  with  some  form  of  boundary  layer 
control.  The  boundary  layer  velocity  profile  on  the  body  is  assumed 
known  and  is  denoted  by  u^.  Likewise  the  wall  shear  stress  on  the  body 
is  T,  . 

D 


Beginning  with  Equation  (16)  of  Reference  [1]  or  Equation  (24)  of 
Reference  [2] 

4 

ACp  =  l  ACp  ,  (3) 

i=l 


where  ACp^  was  shown  to  be  negligible  and 


ACp^  -  2 


fA  1 
-2. 

2 

z  4tt 

[a 

,  ,  2 

<  J 

AU 

00  ! 

rw 


u  r  dr  , 


(4) 


ACp  2  = 


4tt 


ApU 


2  J  rb  Tb  dx' 

00  l 


(5) 


ACp3  =  - 


47rr 

w 


ApU 


n  |  T  dx'  , 
2  |  w 


(6a) 


°-iw  *  91  fc 


1/5 


4/5 


(6b) 


The  approximation  of  Equation  (6b)  results  from  the  assumption  that  a 
flat  plate-like  turbulent  boundary  layer  develops  on  the  duct  wall  with 
an  origin  at  x  =  0.  This  part  of  the  total  correction  [Equation  (3)1 
needs  no  further  discussion  here  because  it  will  be  unaffected  by  body 
boundary  layer  control. 
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In  the  computational  schemes  of  Hoffman  [8]  and  most  others,  the 
wall  shear  stress  is  expressed  in  non-dimensional  skin  friction  coefficient 
form: 


1/2  p  U 


(7) 


Hence,  given  Cfb,  Equation  (5) 
integration,  i.e.. 


may  be  computed  directly  by  numerical 


ACp2 


dx' 


(8) 


The  pressure  coefficient  correction  ACp^  is  evaluated  for  a  body 
with  boundary  layer  control  by  first  writing 


r 


I  = 


u  rdr  = 


r  +6  r  -6 

b/-  b  j  2  w 

ub  rdr  +  “v 


r 

rw 


rdr  + 


2  . 

u  rdr 
w 


(9) 


rb+6b 


r  -6 
w  w 


The  assumption  of  flat  plate  turbulent  boundary  layer  development  on  the 
duct  wall,  reduces  Equation  (9)  to: 


1  =  2  UV2  [(rw  -  V2  "  (rb  - 


2  7  2 

6.)  n  +  i  Uy^r  6 
b  9  v  w  w 


7_ 

16 


(uvV2  +  uv2 


2 


rdr 


(10) 


where  u^/u^  would  be  the  non-dimensional  velocity  profile  for  the  body 
with  (or  without)  boundary  layer  control.  It  is  assumed  here  that 
uv  =  u  (edge  velocity)  and  is  negligibly  affected  by  the  tunnel  wall 
meaning  that  the  velocity  profiles  may  be  calculated  using  Cpj.  Also, 

6  is  the  distance  from  the  body  surface  to  where  u  =  0.99  uy.  Elimination 
of  uv2  from  Equation  (10)  requires  the  continuity  equation: 


r 

A  U  =  2tt 
o  ® 


j*  urdr 


(11) 


t 


t 


# 


i 
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Writing 


A  U  /  2u 
o  00 


'  b 

w  r  w 

u,  rdr  +  u 

i  rdr  + 

b 

V  j 

J 

u  rdr 
w 


V5b 


r  -6 
w  w 


(12) 


one  can  easily  verify  that 


A  U  /2tt  =  ^  u  [  (r  -  6  )2  -  (rL  +  6V)2] 
o  2  v  w  w  b  b  J 


+  u  r  6  -  ~  6  2]  +  u  , 

Vl8ww  15  w  V  J 


r+6 

b  (  b  ru^\ 


V 


rdr 


(13) 


Solving  Equation  (13)  for  u^,  squaring  and  substituting  into  Equations  (10) 
and  (4)  one  finds  that 

-  4+5, 

rt  o  i  /  -J  2  u  f  D 

I  ( r  -6  )  -  (r,+6,  )  +  ^  r  5  --£-6  +2 

1  v  w  w  bb  9  ww  8w 


.  1  2A  ‘ 

rA  i  2  o 


ACpi  =  2 


(uju)  rdr 


b'  V 


7TA 


k  2  .  7 


14  .  2  .  „  rb^5b 


(r  -6  )  -  (r,+6,  )  +tt6  -~6  +2 

w  w  b  b  4  w  w  15  w 


(ub/uv)  rdr 


Equations  (8)  and  (14)  are  seen  to  require  numerical  integration  of  the  skin 
friction  coefficient  and  boundary  layer  velocity  profile,  respectively. 


EXAMPLE 


To  illustrate  the  effect  of  one  type  of  body  boundary  layer  control  on 
the  pressure  distribution  of  an  axisymmetric  body  which  operates  in  the  closed 
cylindrical  working  section  of  a  water  tunnel,  some  numerical  calculations 
based  on  the  preceding  analysis  will  be  presented.  The  body  selected  for 
this  example  is  one  designed  for  laminar  flow  control  through  distributed 
surface  suction.  The  body  is  2.900  m  long  and  has  a  cylindircal  mid  section 
0.318  m  in  diameter.  The  calculations  to  follow  assume  it  will  operate  in 
a  water  tunnel  test  section  which  is  1.22  m  in  diameter.  The  nose  of  the 
body  was  designed  by  Smith  [10]  and  is  described  by  equations  of  the  Lemniscate 
type  [11].  In  particular,  the  nose  contour,  as  defined  back  to  the  tangent 
point  where  the  cylindrical  mid  body  begins  is  given  by: 


x  =  3.6  D  (1  -  1.241  /sin  20  cos6)/2  (15a) 

max 

r,  ■  D  (1  -  2.149  /sin  20  sinG)/2 
b  max 


.(14) 


(15b) 
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where  0  _<  0  <_  30°.  The  tangent  point  defines  the  nose  length  and  is  at 
x  =  0.572  m.  The  tailcone  originates  at  x  =  2.134  m  and  is  described  by 
a  6th  order  polynomial.  Because  the  subject  calculations  are  not  accurate 
over  the  afterbody  [1,2]  region  there  is  no  need  to  give  the  mathematical 
details  of  the  afterbody  here. 

Researchers  at  the  Naval  Ocean  Systems  Center  in  San  Diego,  CA  have 
performed  boundary  layer  stability  analysis  (unpublished)  on  this  particular 
body  design  and  have  concluded  that  the  suction  distribution  given  in 
Figure  2  should  result  in  full  laminar  flow  over  the  body  for  an  extensive 
speed  range.  Here,  the  suction  coefficient  Cq  is  defined  in  the  usual 
manner  as  un/ue,  where  un  is  the  fluid  velocity  normal  to  the  body. 

Typical  laminar  boundary  layer  velocity  profiles  are  given  for  two 
velocities  and  two  different  axial  locations  in  Figures  3  through  6. 

These  profiles  were  computed  using  Hoffman's  computer  code  [8]  which  uses 
the  potential  flow  in-tunnel  body  pressure  distribution  calculated  from 
the  program  of  Friedman  [4] . 

The  calculated  potential  flow  and  viscous  flow  (with  suction)  pressure 
distributions  for  this  body  are  given  for  three  velocities  in  Figures  7 
through  9.  The  pressure  coefficient  corrections  summarized  by  Equation  (3) 
can  be  seen  on  these  figures  as  the  difference  between  the  dash-dot  line 
and  the  dashed  line.  As  expected,  the  correction  is  most  important  at 
low  speeds.  The  effect  of  surface  suction  on  these  viscous  flow  pressure 
distributions  is  demonstrated  in  Figures  10  through  12.  The  "without 
suction"  results  are  based  on  laminar  boundary  layer  velocity  profiles 
computed  with  Cq  =  0.  It  must  be  emphasized,  however,  that  at  these 
velocities  and  lengths,  the  boundary  layer  will  almost  certainly  be  turbulent 
over  most  of  the  body  length  with  no  surface  suction.  Since  turbulent 
boundary  layers  are  thicker  than  laminar  ones,  the  differences  between  the 
pressure  distributions  of  Figures  10  through  12  are  considered  the  minimum 
to  be  expected. 


CONCLUSIONS 


It  has  been  the  purpose  of  this  note  to  extend  an  analysis  [1,2]  for 
the  effects  of  viscous  flow  on  in-duct  body  pressure  distributions  to 
account  for  the  non-conventional  boundary  layer  development  over  bodies 
with  some  form  of  boundary  layer  control  such  as  heat,  suction,  blowing, 
or  cooling.  The  resulting  formulation  is  similar  to  that  for  natural 
flat  plate-like  body  boundary  layer  development  but  is  more  complicated 
because  the  boundary  layer  velocity  profiles  for  the  body  with  boundary 
layer  control  need  to  be  specified  at  all  axial  stations  of  interest  and 
then  numerically  integrated.  Aside  from  this  numerical  inconvenience, 
the  computations  are  straightforward.  An  example  was  given  for  a  typical 
body  of  revolution  with  surface  suction.  Based  on  the  computed  pressure 
distributions  for  this  body  operating  in  a  1.22  m  diameter  duct,  the 
differences  between  the  pressure  coefficients  with  and  without  suction 
were  seen  to  be  significant;  typically,  on  the  order  of  10  percent  at 
Cp  a  -0.2.  Although  examples  were  not  presented  for  other  types  of 
boundary  layer  control,  similar  differences  would  be  expected. 


-11- 


10  August  1981 
GCL: cag 


The  usefulness  of  the  current  analysis  appears  to  be  in  the  area  of 
laminar  flow  control  experiments,  where  a  given  system  might  be  evaluated 
in  a  wind  or  water  tunnel.  Because  of  the  complexity  of  such  bodies  it 
would  prove  difficult  to  instrument  them  with  static  pressure  orifices 
to  obtain  precise  pressure  distribution  data.  (Precise  data  are  required 
because  transition  locations  are  very  sensitive  and  known  to  move  with 
only  slight  changes  in  pressure  gradient.)  Since  these  data  are  required 
in  order  to  execute  properly  most  transition  prediction  methods,  a  method 
such  as  this  one  would  seem  to  be  required  in  order  to  account  for  some 
of  the  many  variables  that  enter  into  such  an  experiment.  Unfortunately, 
there  are  no  currently  available  experimental  pressure  distributions  for 
bodies  with  boundary  layer  control,  so  verification  of  the  current  analytical 
extension  of  the  theory  of  References  [1  and  2]  cannot  be  performed  at  this 
time.  However,  it  is  emphasized  that  the  work  documented  in  References  [1 
and  2]  does  include  experimental  results  for  a  body  without  boundary  layer 
control  that  do  indeed  verify  the  theoretical  approach  taken  and  assumptions 
made. 
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Figure  1 
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a)  Potential  Flow  Over  a  Body  in  an  Infinitely  Long 
Cylindrical  Duct;  b)  A  Schematic  Description  of  Boundary 
Layer  Development  on  the  Body  and  Duct  Wall,  Noting  That 
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Figure  2.  Theoretically  Established  Suction  Distribution  for  Laminar 
Flow  Control  of  Body  Considered  in  the  Numerical  Example 
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Figure  5 


Body  Boundary  Layer_Velocity  Profiles  Computed  With  and 
Without  Suction  at  x  *  1.768  and  6.1  m/sec 
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Figure  6 
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Figure  7.  Computed  Potential  Flow  Pressure  Distributions  of  Body 
Operating  in  Free  Stream  and  in  a  1.22  m  Diameter  Duct. 
Also,  the  In-Duct  Result  Corrected  for  Viscous  Effects  When 
Suction  is  Applied  and  for  -  6.1  m/sec 
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Figure  8.  Computed  Potential  Flow  Pressure  Distributions  of  Body 
Operating  in  Free  Stream  and  in  a  1.22  i  Diameter  Duct. 
Also,  the  In-Duct  Result  Corrected  for  Viscous  Effects  When 
Suction  is  Applied  and  for  U  =  12.19  m/sec 


10  August  1981 
GCL: cag 


V 


Figure  10.  Comparison  of  the  In-Tunnel  Viscous  Flow  Pressure 

Distributions  With  and  Without  Body  Suction  at  6.1  m/sec 


-23- 


10  August  1981 
GCL: cag 


BODY  LENGTH  (m) 


Figure  11.  Comparison  of  the  In-Tunnel  Viscous  Flow  Pressure 

Distributions  With  and  Without  Body  Suction  at  12.19  m/sec 
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